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Abstract. We investigate ground-state properties of two correlated-hopping electron 
models, the Hirsch and the Bariev model. Both models are of recent interest in the 
context of hole superconductivity. Applying the Lanczos technique to small clusters, we 
numerically determine the binding energy, the spin gaps, correlation functions, and other 
^ . properties for various values of the bond-charge interaction parameter. Our results for 

small systems indicate that pairing is favoured in a certain parameter range. However, in 
contrast to the Bariev model, superconducting correlations are suppressed in the Hirsch 
^> ' model, for a bond-charge repulsion larger than a critical value. 
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Since the discovery of high-T c -superconductivity strong efforts have been made 
to explain the pairing mechanism with the aid of microscopic models for strongly 
correlated fermions, the Hubbard and t-J model being the most popular among 
them. In one dimension, exact results can be derived for the Hubbard model 
and for the supersymmetric t-J model by applying the Bethe ansatz |1J]. The 
main interest is the question whether and how collective behaviour can lead to 
a compensation of the Coulomb repulsion between the electrons. As a possible 
explanation, models with correlated-hopping interactions are a subject of current 
c3 ' research &4T 



In 1989, Hirsch proposed a model for the description of oxide superconduc- 
tors by considering the holes in a nearly filled band as the charge carriers. The 
Hamiltonian contains among other contributions a correlated hopping interaction, 
i. e. a bond-charge repulsion. A modified version of Hirsch's model, the Bariev 
model, has been solved by Bethe ansatz [|3] . Although it only contains one half of 
the bond-charge interaction terms of the original model, it is expected to maintain 
its basic qualities because it also takes into account the modification of the hop- 
ping amplitude by the presence of a particle with opposite spin. Considered as an 
electron model, the bond-charge repulsion leads to the formation of Cooper pairs 
of holes, a process which is favoured at small hole density ||] . 

In this paper we compare exact diagonalization results for the two models ob- 
tained by using the Lanczos technique Jl9| . Although the systems treated are small 



(up to 16 sites) and despite the finite size effects which render the interpretation of 



1 



-15 







1 








- 

*'./ / : 
,* /* : 




\ 


"x- 


■*• 




..* 


* : 








X 









5 10 15 20 

N 

Figure 1: Ground-state energy, L — 12 sites, of the Hirsch (+) and the Bariev 
model (x) with PBC. At = —1,-0.2,0.2,1 from bottom to top. The analytic 
result for At = 1 is indicated as full line. 



the results difficult, comparison with some exact results for the Bariev model (in 
the thermodynamic limit), obtained from the Bethe ansatz solution, shows good 
agreement. 

In the following we shall consider a one-dimensional system consisting of N = 
iVf + iVj_ < 2L itinerant electrons on a closed chain of L sites with periodic (PBC) 
or antiperiodic (ABC) boundary conditions. Electron hopping is possible between 
nearest neighbour sites, but the hopping matrix element is modified by the presence 
of electrons with the opposite spin direction on the sites involved in the hopping 
process. 

The Hamiltonian of the Hirsch model reads 

Hn(At) = -J2J2 (4,^+^ + c t+i,° c i>°) ( 1 ~ -y ( n ^- CT + ) > (!) 

j = l cr=±l ^ ' 

and the Bariev model is described by the Hamiltonian 

L 

H B (At) = -J2 E ( c t°j+^ + '•;•!.'.'>) ( x " ten i+ 2$L,-) ' ( 2 ) 

3=1 <r=±l 

where c^ CT and c JiCr , respectively, denote the creation and annihilation operators of 
an up-electron (a = 1) or down-electron (a — — 1) at site j. rij y(T (which equals 
or 1) is the number of particles with spin a at site j. As can easily be seen, the 
interaction parameter At is positive in the repulsive model because it leads to 
a decrease of the hopping probability (for At < 2), if an electron with opposite 
spin direction is present. In this sense there is an additional correlation between 
the up-spin and down-spin hopping processes. 

We start our exact diagonalization study with the discussion of the ground- 
state energy. The energy difference between the Hirsch and the Bariev model is 
very small for small At. As shown in fig. |I|, the energy curves, as functions of 
the particle number N, coincide for At = —0.2 and At = 0.2. They are also 
in good agreement for At = —1, but the curves spread apart for At = 1, where 
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Figure 2: Spin gap of the Hirsch model (left) and Bariev model (right) for n = 3/2 
and L = 4, 8, 12, 16 from top to bottom. 



the difference between the bond-charge interaction terms in the two models be- 
comes relevant. The figure also indicates that the energy difference is small at low 
electron and low hole density for all values of the interaction parameter At. For 
comparison, we included the exact result H] for the Bariev model at At = 1 as a 
full line, which is = — ~(1 + |)sin(^W). The exact diagonalization results for 
L = 12 all lie on this curve. 



In the following, we discuss properties of the Hirsch and Bariev model for small 
hole filling. In the case of quarter filling of holes (i. e. electron filling n = 3/2), we 
can compare results for systems with L = 4, 8, 12, and 16 sites. 
The spin excitation energy is defined as A s = E(S Z = 1) — E(S Z = 0), where 
E(S Z ) is the ground-state energy in the subspace with fixed S z ; it corresponds 
to the lower critical magnetic field h c which is necessary to create a non-zero 
magnetization in the system. In the absence of a spin excitation gap, h c vanishes, 
as is the case in the repulsive Hubbard model and the supersymmetric t-J model. 
On the contrary, the attractive Hubbard model and the Bariev model have a gap 
in the spin excitations (see |5], |5], |2D| and references therein). 



In fig. g we show the results for A s for both models. We choose PBC for L = 16 
and 8 and ABC for L =12 and 4, because the results facilitate the extrapolation 
into the thermodynamic limit if the Fermi wavenumber hp = nir/2 of the free 
infinite system is one of the possible k numbers of the finite system |21[|. We 



verified that there is already good agreement between the exact diagonalization 
results for L = 12 and the analytic solution A s = | sin — 2+^ cos (1+^) 

for the Bariev model at At = 1. The results strongly suggest that the spin gap is 
maximal at At = 1 in both models and that the spin gap is vanishing in the Hirsch 
model for At > 1.5 in contrast to the Bariev model where it remains non-zero. The 
vanishing of A s at At = 2 in the Hirsch model is evident for any system size and 
boundary condition and a large range of densities. It is due to the high degeneracy 
at that parameter value where the Hirsch model can be treated analytically and 



the ground-states are of the 7/-pairing type M, ITTJ . Moreover, results for other hole 



densities indicate that the vanishing of the spin gap in the Hirsch model is not 
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Figure 3: Spatial singlet pair correlations Go (°), G\ (□), and triplet pair correla- 
tions G2 (a) of the Hirsch model at At = 1 (left) and At = 1.7 (right) for L = 16, 
n = 3/2. 
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Figure 4: Same as in fig. [3] for the Bariev model. 



restricted to the band filling of n — 3/2. 

Whenever there is a gap in the spin excitation spectrum, spin correlations decay 
exponentially so that other correlations, such as pair correlations are more likely 
to dominate. We calculate the singlet pair correlation functions for on-site and 
next-neighbour pairs and the triplet pair correlation function, defined as 

G (r) = (P+(r)P (0)>, 
d(r) = (P+(r)P s (0)>, 
G a (r) = (P 4 + (r)P(0)), 



with the pairing operators 



P+(r) 



P s + (r) 



M L r| c r+1,TJ' 

1 

V2 V 



-P* + ( r ) — /^( C ^T C ^+l,i + C rl C r+l,l)- 



When the number of doubly occupied sites is reduced in the ground-state of 
the interacting system — as compared to free fermions — Gi(r) is expected to be 
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Figure 5: Fourier transform G(k) of the pair correlation functions of the Hirsch 
model (left) and Bariev model (right), at k — 0. G (o), G\ (□), G 2 (a). 



larger than Go(r), although both functions should be equivalent with respect to 
their decay in the infinite system. We compare results for both models for At = 1 
and At = 1.7. Fig. |3| shows that at At = 1 both singlet correlation functions of 
the Hirsch model decay slowly to a non-zero value (which nevertheless can be zero 
in the infinite system, when finite size effects are absent). The triplet correlations 
decay very rapidly to zero. On the other hand, all pair correlations decay fast to 
zero at At = 1.7 where there is no spin excitation gap. This is different for the 
Bariev model where the decay is slower (fig. |j). We mention that the spin operator 
S 2 does not commute with the Hamiltonian Tig (0) so that S is no good quantum 
number in the Bariev model, in contrast to the Hirsch model. 

To gain some insight into the dependence of the pair correlation functions on 
the bond-charge interaction parameter At, we also show their Fourier transform 
Gj{k) = Y2 r cos(fcr) (j = 0, 1, 2) for k = which is simply the sum of the 
correlations over all spatial distances. Although it is not possible to discuss the 
decay by means of this function only, it helps to determine at which values of the 
interaction pair correlations are favoured. Fig. [| indicates that singlet pairing is 
strong at At ~ 1 in both models, i. e. where the spin gap is large. This tendency 
is diminished drastically in the Hirsch model at parameters At > 1.1. On the 
other hand, these functions are decreasing more slowly in the Bariev model. The 
curves of Go and Gi have a similar shape and only differ by magnitude. Triplet 
correlations are small in both models and vary little as functions of At. 

To complete the comparison of the two models, we calculate the binding energy 
E B (N) = E(N + 2) + E(N) - 2E(N + 1). To suppress the finite size effects, we 
choose twisted boundary conditions (TBC), which means that the wavefunction 
at site and site L is related by ^(L) = e i</3 \I/(0), and take the average of Eb 
over various values of the twisting angle <p. To estimate the finite size effects, we 
also show the maximum/minimum Eb obtained in this way (fig. |GJ). The binding 
energy is negative in the Hirsch model for At < 1.4, but becomes positive for 
larger values of At. In the Bariev model, the averaged Eb is zero at At = and 
At = 2 and negative for all other At, with minimum at At = 1. For At > 1.7, the 
finite size effects in the Bariev model are too large to decide if binding occurs or 
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Figure 6: Binding energy of the Hirsch model (left) and Bariev model (right), 
n = 20/12, averaged over TBC (full line). The broken lines estimate the finite size 
deviations. 



not, but we have strong evidence that binding is suppressed in the Hirsch model 
for At > 1.5. 

We conclude that our calculation of various ground-state properties of finite 
clusters at small hole doping indicates that the Bariev model and the Hirsch model 
show similar behaviour as long as the bond-charge repulsion At is small. Pair 
correlations are strongest in both models for At ~ 1. For At > 1.5, there is a 
change in the ground-state properties of the Hirsch model, and superconducting 
correlations are strongly suppressed for large At, which is in accordance with 
mean field results In the Bariev model, there is no such evident change. 
We finally propose a simple argument which might help to understand why the 
parameter of At ~ 1.5 can be a phase boundary in the Hirsch model: The analytic 
diagonalization of the local interaction (acting only on two neighbouring sites) of 
d) shows that (for 1 < At < 2) = \a2) + \2a) and |tf 2 ) = |||) - \iX) + 
1 20) + 1 02) (0, a, 2 denoting the empty, singly, and doubly occupied site) are 
local eigenvectors with eigenvalues Ei — 1 — At and E 2 = — 2 + At ||12|| . Since 
Ei = E 2 for At = 3/2, there is a level-crossing in the local ground-state at that 
parameter, whereas this effect is absent in the Bariev model. This criterion might 
be of importance for the global ground-state, too. 
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